Effect of non-zero constant vorticity on the nonlinear resonances 

of capillary water waves 
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The influence of an underlying current on 3- wave interactions of capillary water waves is studied. 
The fact that in irrotational flow resonant 3-wave interactions are not possible can be invalidated 
by the presence of an underlying current of constant non-zero vorticity. We show that: 1) wave 
trains in flows with constant non-zero vorticity are possible only for two-dimensional flows; 2) only 
positive constant vorticities can trigger the appearance of three-wave resonances; 3) the number of 
positive constant vorticities which do trigger a resonance is countable; 4) the magnitude of a positive 
constant vorticity triggering a resonance can not be too small. 
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1. Introduction. In this Letter we investigate the 
effect of a current on the dynamics of nonlinear capillary 
water waves, a problem of fundamental importance 
due to the ubiquity of currents at sea [8]. The most 
common force for creating water waves is the wind and 
wind generated capillary waves play a prominent role 
in the development of waves on water surfaces that are 
flat in the absence of wind. Indeed, capillary waves 
generate surface roughness allowing a better grip of 
the wind. This leads to the subsequent development of 
capillary-gravity and gravity waves as with increasing 
wave amplitude gravity becomes the dominant restoring 
force replacing surface tension [15, 18, 20]. In coastal 
navigation the important question arises whether the 
presence of an underlying current can be detected by 
investigating solely phenomena at the water's surface. 
Vorticity is adequate for the specification of a cur- 
rent. A uniform current is described by zero vorticity 
(irrotational flow), while the simplest example of a 
non-uniform current is that of tidal flows, which can 
be realistically modeled as two-dimensional flows with 
constant non-zero vorticity, with the sign of the vorticity 
distinguishing between ebb/tide [7]. Notice that in 
linear systems waves of different frequencies do not 
interact due to the superposition principle, while in 
a nonlinear system the lowest-order nonlinear effect 
(with respect to an expansion in the wave amplitudes) 
is the resonant interaction of three waves of different 
wave-vectors and frequencies. Resonant interactions 
can profoundly affect the evolution of waves by making 
significant energy transfer possible among the dominant 
wave trains, accounting thus for wave patterns that 
are higher and steeper than linear wave theory would 
predict and providing insight into the effects of weak 
turbulence [19]. It is known that capillary waves in 
irrotational flow do not exhibit exact 3-wave resonance 
[11]. We will show that the presence of an underlying 
current of constant vorticity can lead to the excitation 
of nonlinear resonances, but only for special vorticities. 



2. Dispersion function. Our purpose is to show that 
currents strongly affect the dynamics already at the level 
of capillary waves. To emphasize the wave-current inter- 
actions we consider a setting with a flat bed in order to 
minimize the effect of the shape of the shoreline and of 
the water bed on such flows. In this context we investi- 
gate whether or not three-wave resonances are possible 
among capillary waves. The dispersion relation for exact 
solutions to the governing equations for capillary-gravity 
water waves propagating at the free surface of water of 
constant density p = 1 with a flat bed and in a flow with 
constant vorticity Q is 
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where d > is the average water depth, g is the gravi- 
tational constant of acceleration and a > is the coef- 
ficient of surface tension, and the notation k is used for 
the modulus of the wave vector, k = |k|. The case of 
capillary waves corresponds to g — ► and that of grav- 
ity waves to o — ► cf. [4, 22, 23]. Moreover, capillary 
waves being short waves, it is appropriate to take the 
short wave limit kd — > oo of the above dispersion rela- 
tion; as tanh(fcd) — > 1 we obtain the dispersion relation 
for small-amplitude capillary waves as 
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In the case of irrotational flows (with zero vorticity), the 
dispersion relation reads [23] 
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As it was shown in [11], exact three- wave resonances are 
prohibited for (3) by purely kinematic considerations. 
Below we show that in the case of non-zero constant 
vorticity, resonances are possible but only for a countable 
number of positive vorticities. 
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3. Magnitudes of vorticity. First, we show that the 
resonance conditions 

w(ki) + w(k 2 ) =w(k 3 ), ki+k 2 =k 3 (4) 

do have solutions for the dispersion function (2). Let us 
introduce the variable £ = Q/2^/a and rewrite (4) as 

Vc 2 + fc? + ^C 2 + fci = ^C 2 + fc| - C, k 1 +k 2 = k 3 . 

Our primary goal is investigating whether for a given 
current strength f2 three- wave resonances are possible for 
wave vectors with integer coordinates. Since capillary 
wave trains propagating in flows of constant vorticity 
are only possible if the flow is two-dimensional (see Sec- 
tion 5), only the case of one-dimensional wave- vectors, 
dimk = 1, is relevant. Now k = k and (4) reads 



Vc 2 + k\ + ^Je + ki = ^Je + h - c, (5) 

h+k 2 = k 3 , (6) 

and we will study the magnitudes of vorticities corre- 
sponding to integer wave- vectors fci, fc 2 , fc 3 > 0. To in- 
vestigate (5), (6) it is appropriate to view it as an equation 
in the unknown £ ^ 0, with k±, fc 2 , fc 3 > given integers. 
Let us first simplify notation by denoting 

k\ = a, fcf = /3, fc|= 7 , E = a + p~ 1 . (7) 

Since irrespective of the sign of £ both sides of (5) are 
positive, squaring them we see that (5) is equivalent to 



E + 2C V /C 2 +^ = -2v/(C 2 + a)(C 2 + J)- (8) 
Once we obtain a formula for £ 2 , equation (8) yields 

£ + 2 v /(C 2 +a)(C 2 +/3) 



c = - 



But (8) squared yields 
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(9) 



E 2 - 4( 2 E - 4a(3 = -4#C\/C 2 + 7, 
which, squared again, leads to 

(E 2 - 4a/?) 2 = 

8E C {2E 1 + E 2 - 4a/?) = 8E ( 2 [(a - /?) 2 - 7 2 ] , 

with the last equality a simple consequence of the defini- 
tion of E. Thus 



(E 2 - Aa/3) 2 



1 (a + /3- 7 )[(a-/3) 2 ~7 2 1' 



(10) 



Obviously, (10) can be regarded as a 2-parameter series 
of solutions for (5)- (6). Indeed, for any wave- vectors 
fci , fc 2 there are at most two vorticities VL generating a 
three- wave resonance of the form (5)- (6). For instance, 



h = k 2 = 1 yield fc 3 = 2, a = /? = 1, 7 = 8, E = 
—6, ( 2 = 1/3 and the corresponding vorticity ft has to 
satisfy ft 2 = 4ct/3. Now we have to gain more under- 
standing about the possible sign of the vorticity. To for- 
mulate our theorem on the sign of vorticity, let us first 
introduce a new variable 



(11) 



Theorem 1 (on the sign on vorticity). Three wave- 
vectors fci, fc 2 , fc 3 > 1 are solutions of (5), (6) if and only 
if 5 > 0, in which case the vorticity f2 = 2 v / cr£ is positive, 
with £ given by (9). 

Proof. Two observations can be made immediately. 
Firstly, as it was shown in [11], the equality 5 — can be 
transformed to the particular cubic case of Fermat's Last 
Theorem and therefore has no integer solutions. Sec- 
ondly condition (6) yields 



kl > k\ 



k 2 



(12) 



and therefore E < 0. 

We first claim that if fci, k 2 , fc 3 > 1 are solutions of 
(5), (6), then 5 > 0. Indeed, notice that the inequality 



3fc 2 fc 2 + 3fcifc 2 = 3fcifc 2 (fci + k 2 ) > 6kik 2 y / hk 2 



yields 3fc 2 fc 2 + ik Y k 2 2 > 2^Jk\k\. If (6) holds, adding 
k\ + k 2 to both sides of the above inequality we get 
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that is, S > 0. 

Let us now prove that for S > we always have a 
solution C > and no solution £ < 0. The statement 
about the positive solution follows at once by noticing 
that the function C >-> VC 2 + a+ y/( 2 +/3- VC 2 +7 + C 
has limit +00 for y — > 00, being strictly increasing on 
(0, 00) since its derivative satisfies 
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1 > 1, 



for ( > 0, as 7 > /? in view of (12). Therefore this 
function has a zero on (0, 00) if and only if its value at 
( = is strictly negative, that is, if and only if S > 0. 
Now notice that we know by (9) and (10) that if a solution 
exists, it is unique. This concludes the proof. 

Corollary 1 The set of all resonant triads covered by 
(5), (6), is generated by a countable number of positive 
vorticities. 

More precisely, given the wave- vectors fci, k 2 > 1 we 
define fc 3 via (6), and compute 



E = 
E 2 



-3fcifc 2 (fci + k 2 ), 
4a/? = k\k\ (9k 2 



9fc 2 + 14fcifc 2 ), 
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finding also that -f 2 — (a — (i) 2 equals 

kik 2 (6fc 4 + I5kfk 2 + 22k\k\ + 15/cifcf + 6k%). 

Knowing that £ > 0, we infer from (10) that a three-wave 
resonance occurs only if f2 is given explicitly by 

^kl ,2 kl /2 (9k 2 + 9k% + 14Ma) 

v/6(fci + k 2 )(6kf + 15k'(k 2 + 22k\k 2 + 15fc^| + 6fc 4 ) ' 

(13) 

For fci > k 2 we see that the denominator is bounded from 
above by 16^/3 k^ 2 and from below by I6y/3k^ 2 . Since 

9k\ + 9k\ + Ukik 2 = 9(h - k 2 ) 2 + 32fcifc 2 > 32fcifc 2 , 

a lower bound for the numerator is 32 ^/a k^ 2 k^ 2 , with 
the evident upper bound 32 y/a k^ 2 k^ 2 ■ Therefore 

V 3 1 k 2 ~ ~ V 3 2 

Corollary 2 Three-wave resonant interactions do not 
occur in flows with sufficiently small constant vorticity, 
the minimal magnitude of a resonance generating vor- 
ticity being Q, rn i n — 2V^/V3. 

4. Structure of resonances. In the previous section 
we have shown that for any two wave- vectors ki,k 2 , a 
corresponding resonant triad (fci, k 2 ,ki + k 2 ) will be gen- 
erated by the vorticity given by (13). The inverse prob- 
lem — to compute all resonant triads generated by a 
given vorticity — is much more complicated. Already 
for one fixed wave-vector, this problem is equivalent to 
finding rational points on an elliptic curve. On the other 
hand, it is relatively easy to compute numerically, in a 
some finite spectral domain, the set of all resonant triads 
and the corresponding magnitudes of the vorticity. An- 
other reasonable question is whether different resonant 
triads can be generated by almost the same magnitude 
of vorticity, with some accuracy e because these type of 
resonances might enrich the cluster structure. 

The structure of resonance clusters has been investi- 
gated numerically, in the spectral domain k\,k 2 < 100. 
All exact solutions for £ 2 , that is, solutions with e = 
have one of two forms: 

1. Each resonance triad of the form (fci, fci, 2k\) with 
ki = 1,2, ...,100 can be generated by only one vortic- 
ity; two different triads of this form, (fci,fci,2fci) and 
(k 2 ,k 2 ,2k 2 ) with ki ^ k 2 can only be generated by two 
different vorticities. 

2. Each resonant triad of the general form (fci 7 k 2 ,ki + 
k 2 ) with k\ ^ k 2 has the corresponding symmetrical res- 
onant triad (fc2,fci,fci + k 2 ). Both symmetrical triads 
are generated by the same vorticity. This fact can be 
observed immediately from the form of (13) which is in- 
variant under the transformation [fci i— » k 2l k 2 ^> fci]. 



It appears that a given vorticity can generate not more 
than two exact resonant triads. The situation does not 
change substantially when e ^ 0: again, not more than 
two resonant triads appear to be generated by almost the 
same magnitude of vorticity. For example, 



e = 10" 1 => (1,9,10), 


c 2 


= 4.824; 


e = 10" 1 => (2,3,5), 


c 2 


= 4.764; 


e = 10~ 2 => (2,14,16), 


c 2 


= 29.5622; 


e = 10~ 2 => (4,5,9), 


c 2 


= 29.5612; 


e = 10~ 4 => (2,88,90), 


c 2 


= 196.212144; 


e = 10~ 4 => (3,40,43), 


c 2 


= 196.212121; 


The only difference with the case e 


= is that now the 



two co-existing triads are not symmetrical anymore. In 
Table I the distribution of the resonance generating vor- 
ticities is given. One can see immediately that more that 
30% of all vorticities are relatively small, in the sense 
that they satisfy the condition = 2^ ^fa < 2 ■ lO 2 -^- 
(first row on the left in the Table I). The number of 
resonance generating vorticities in the partial domains is 
decreasing exponentially with the growth of the general 
spectral domain. Interestingly enough, new asymmetri- 
cal triads appear for e of the order of 10~ 4 or bigger, i.e. 
e is bounded from below (cf. detuned frequency balance 
in the irrotational case [12, 21]). 



Magnitude of £ 


Number of O 


Magnitude of £ 


Number of O 


C< i 


3477 


4 < C < 5 


494 


1 < C < 2 


1146 


5 < C < 6 


407 


2 < C < 3 


781 


6 < C < 7 


363 


3 < C < 4 


609 


7<C 


2723 



TABLE I: The distribution of the magnitudes of the reso- 
nance generating vorticities is presented, in the spectral do- 
main ki, fe < 100. The notation ( = ( 2 ■ 10 4 is used. 



Generically, to describe the dynamics of the resonances 
in the rotational case, two steps have to be performed. 
Firstly, nonlinear evolution equations governing the 
rotational case should be derived. Afterwards the 
standard procedure would be to construct all resonant 
clusters met in the solution set above and write out 
explicitly the corresponding dynamical systems [13]. 
Secondly, the interaction coefficient has to be computed 
which, of course, will be different from the coefficient 
known for irrotational flows. Computing the explicit 
form of the interaction coefficient for a given wave 
system is a non-trivial technical problem demanding 
tedious computations which can not be automatized in 
the present state of symbolic programming [14]. 

5. Flows with constant vorticity. Let us recall the 
governing equations for capillary water waves propagat- 
ing at the free surface of a layer of water above a flat 
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bed [9]. In the fluid domain of average depth d > 
bounded above by the free surface z = d + r/(x, y, t) and 
below by the fiat bed z = 0, the velocity field u(x,i) and 
the pressure function P(x, t), where x = (x,y,z) and 
u = (u\,U2,Us), satisfy the Euler equation 

u t + (u- V)u + VP = 0, (14) 

as well as the equation of mass conservation 

V-u = 0, (15) 

expressing homogeneity (with constant water density 
p = 1). These equations are coupled with the bound- 
ary conditions 

P = aH, u 3 = r)t + uiT] x +U2r]y, (16) 

on the free surface z = d + f](x, y, t), the constant a > 
being the surface tension coefficient and H being twice 
the mean curvature, 



H(x,y,t) = 



(! + vDVyy - ^VxVyVxy + (1 + vl)Vx 



On the flat bed z = we require 
u 3 = 0. 



(17) 



The governing equations for capillary waves are (14)-(17). 
These equations are well-posed [6]: starting with an ini- 
tial surface profile t)(x, y, 0) with Sobolev regularity H 6 5 
and an initial velocity field u(x, 0) with Sobolev regular- 
ity H 5 - 5 satisfying (15), for some time t € [0,T) with 
T > there exists a unique solution (u, r/, P) of the 
system (14)-(17). In our context it is important to keep 
track of the vorticity fi(x, t), obtained at any time t from 
the velocity profile as V x u = fi. With this purpose, 
notice that using (14)-(17) one can derive the vorticity 
equation [16] 



fit + (u • V) fi = (fi • V) u. 



(18) 



To investigate further the vorticity, it is useful to intro- 
duce the flow map x $(x, t): this map advances each 
particle in the water region from its position x at time 
t = to its position $(x, t) at time t. For fixed t, <E> is 
an invertible smooth mapping and from (14), (15), (18) 
one can infer that [16] 



fi($(x, t), t) = J(x, t) fi(x, 0), 



(19) 



where J(x, t) is the Jacobian matrix of the flow map. 
An immediate consequence of (19) is that in three- 
dimensional flows a particle which has no vorticity never 
acquires it and conversely, a particle which is moving ro- 
tationally will continue to do so. 

Theorem 2 (on the dimension of flows with constant 
vorticity). Capillary wave trains can propagate at the 



free surface of a layer of water with a flat bed in a flow 
of constant non-zero vorticity only if the flow is two- 
dimensional. 

Proof. A wave train is a periodic surface wave which 
propagates without change of shape at constant speed 
c > in a fixed direction, say, that of the rr-coordinate, 
and which is unchanged in the y-direction (horizontal and 
orthogonal to the wave propagation direction). Notice 
that for a two-dimensional water flow (in our setting, 
independent of the y-coordinate) , we have fii = fi3 = 
so that (fi • V) u = fi 2 u y = and therefore (18) implies 
that fi t + (u • V) fi = 0: the vorticity of each individual 
water particle is conserved as the particle moves about. 
In particular, if initially the vorticity is constant, it will 
stay so. The existence of wave trains is ensured in this 
setting provided the wave speed is given by the dispersion 
relation (1) cf. [24]. 

Conversely, consider the wave train z = d + i](x — 
cf) propagating in a flow of constant vorticity fi ^ 0. 
Then (fi • V) u = by (18), that is, at every instant the 
vector u is constant in the direction of fi. We first claim 
that for non-flat free surfaces, the direction fi has to be 
horizontal. Indeed, if SI3 7^ 0, this in combination with 
(17) would yield that u 3 = throughout the flow. The 
third component of (14) then forces P z = so that 



P(x, y, z, t) = aH(x — ct) 



(20) 



throughout the flow in view of the first equation in (16). 
Moreover, d z U2 = — fii and d z u\ = fi 2 yield 



m(x,y,z,t) = vi{x,y,t) + fi 2 z, 
u 2 (x, y, z, t) = v 2 (x, y, t) - fi^. 



(21) 



Further, from (15) we infer the existence of a function 
ip(x, y, t) satisfying 



1px = ~V2, 1p y = V 1 . 



(22) 



Writing the first two components of (14) on the flat bed 
z = and on z = e > with e > small enough for this 
horizontal plane to be in the fluid domain, we get 

^2tpxy - filV'w = -^2lpxx + filV'xy = 0. 

As fi3 = d x u2 — dym yields 

(d 2 x + d 2 y )ip = -n 3 , (23) 

and (fi • V) u = yields 

filV'zy + ^2^yy + fi2fi3 = "filV'xx - ^2lpxy - filfi3 = 0, 

unless fii = fi 2 = we infer from these relations that 



i>x 



of +fi^ 



2 , o2> ^ x v ~ (->-' _ o2' wy 



, Ipy 



"fif+fif' 
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Therefore, assuming ft 2 + 0| > 0, we get 
ip(t, x,y) = Ay 2 + Bxy + Cx 2 + a(t) x + b(t) y + kit), 
for some functions a, b, k, where we denoted 

_ ^2^3 r> _ ^1^2^3 _ fl 2 fl 3 

~ ^2(^ + ^2)' ~ Of + r2 2> ° - ~ 2 (ft 2 + f>2)- 
From (14) we now infer that 

^ = -b'(t)-Bb(t)+2Aa(t), P y = a'(t)+2Cb(t)-Ba(t), 

and (20) then yields 

P(x,y,t) = a(x-ct) + (3 

for some constants a, (3. But then P = oH is im- 
possible since the right-hand side has by periodicity in- 
finitely many zeros (at least one between two consecu- 
tive crests, by the mean- value theorem). To rule out the 
possibility that SI3 7^ 0, it remains to consider the case 
Oi = f^2 = 0. If this holds, then the functions u\ and 112 
are independent of z in view of (21) and are harmonic in 
(x,y) by (22)-(23). The second equation in (16) yields 
[ui(x,y,t) — c] r/ x (x — ct) = for all x, y real. But the 
function u\ is a harmonic function of (x, y) and taking 
into account the structure of the level sets of a harmonic 
function (they are curves in the plane unless the function 
is constant cf. [2]), this yields u\ = c since r\ x ^ 0. But 
then the first equation in (14) ensures P x = 0, while (15) 
yields that U2 is independent of y. As Q3 — d x U2 — d y U\ 
we have U2(t, x, y) = Q 3 x+fit) for some function /. Now 
the second equation in (14) yields P y = —f'(t) — cfl3. But 
this can be reconciled with P(x, y, t) = aH(x — ct) only if 
H is a constant and then the surface has to be flat [1, 17]. 

We therefore proved that fl 3 = 0. It remains to show 
that fii = and that P and u are independent of y. 
Since Q\d x Ui + Q 2 d y Ui = for i = 1, 2, 3 and V x u = O 
throughout the flow, we first infer that 

ftiui + VL 2 u 2 = h(t) (24) 

for some function h, since all spatial derivatives of the 
left-hand side are zero. Multiplying the first equation in 
(14) by Oi, the second equation by Q2, and adding up, 
we infer from the above relations that 

n 1 P x + n 2 P y = -h'{t) (25) 

throughout the fluid. Notice that once we show that 
tti = 0, then tt 2 ^ as Q, ^ 0, and (Q ■ V)u = en- 
sures that u does not depend on y. Moreover, in this 
case (24) yields that U2 is simply a function of time and 
consequently the second equation in (14) forces P to be 
independent of y. To complete the proof it therefore suf- 
fices to show that fti = 0. Assuming SI 1 7^ 0, the relation 
d x u 2 —d y u\ — ensures the existence of a function ip(x, t) 
satisfying 

Ul = ip X , U 2 = tfiy. 



From V x u = f! we infer 

d v u 3 = ip yz + Six, d x u 3 = ip xz - n 2 , 

so that, adding if necessary a function depending only on 
(z, t) to ip, we have 

u 3 = ip z + Qiy - fl 2 x. 

Denoting £ = Q 2 x — Qiy, we get 

ui = f x , u 2 = tp y , u 3 = fz~ £. (26) 

and (15) becomes 

<Pxx + fyy + Vzz = 0. (27) 

From (24) and (25) we deduce that <p and P admit the 
decomposition 

ip(x,t) = (fix + vy)h(t) + F(£,z,t), (2& 
P(x,t) = -(jix + vy)h'{t) + G(S,z,t), 

for some functions F, G, where we denoted 
Using (26) and (28), we can express (14) as 

n 2 o s [F t + ^ F | + 1 F 2 _ £ Fz + G ] = 0, 

Hi d 6 [F t + ^ p* + \F 2 - £F Z + G] = 0, 
d z [Ft + F| + \ Ft - iF z + G} = (fi? + fi§) F s . 

Under our working assumption Oi 7^ these relations 
imply F# = 0. But since (27) and (28) yield 

F zz + (f! 2 + Q 2 ) F« = 0, 

we also have F zz — 0. Consequently 

F(t, z, t) = f {t) iz + h{t) e + f 2 (t) z + hit) 

for some functions /j, i = 0, 1,2,3. From (26) and (28) 
we now infer that 

'«i(x,t) = -n 2 [/ (t)« + /i(t)] + MM*). 

< u 2 (x,t) = Oi [hit) z + fiit)] +uh(t), 
^ u 3 ( X) t) = [hit) - 1] in 2 x - Cliy) + hit). 

But then equating the coefficient of y on both sides of 
the second relation in (16), stating that 

u 3 = iui — c) r) x ix — ct) on z = d + 77(2; — ct), (29) 

we first infer that hit) = 1 as Hi 7^ by assumption. 
But then U3(x, t) = hit) and now the right-hand side of 
(29) is periodic in x and the left-hand side is independent 
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of x. Therefore rj x = which contradicts the fact that 
the free surface was not fiat. This concludes the proof. 

Theorem 2 allows us to consider for flows of constant 
vorticity only two-dimensional flows propagating in the 
x-direction, for which the vorticity vector takes the form 
(0,^2,0). For such flows it is therefore natural by abuse 
of notation to identify the vorticity vector with its 
second component. For irrotational flows £1 = so that 
a velocity potential exists and methods from harmonic 
function theory (involving the Dirichlet-Neumann opera- 
tor) can be used to transform the governing equations for 
capillary water waves to a Hamiltonian system expressed 
solely in terms of the free surface rj(x,i) and of the 
restriction of the velocity potential on the free surface 
[25]. The absence of a velocity potential for non-zero 
vorticities complicates the analysis considerably cf. [5]. 

6. Discussion. Our main conclusions can be formulated 
as follows: 

• Flows with constant non-zero vorticity admit capil- 
lary wave trains only if they are two-dimensional. 

• Only positive vorticity can trigger the appearance 
of three-wave resonances. 

• The number of positive vorticities which do trigger 
a resonance is countable. 

• The magnitude of a positive vorticity triggering a 
resonance can not be too small. 

It is remarkable that relatively large positive constant 
vorticities are necessary to observe resonant 3-wave in- 
teractions. The fact that resonances can never occur for 
negative constant vorticities substantiates the belief that 
the vorticity of the flow has a big influence on the dynam- 
ics of the surface water waves. In this context the "frozen 
turbulence" observed in numerical simulations [19] for 
capillary waves is just a manifestation of the fact that 
resonances are absent. 

To understand the nonlinear resonance dynamics 
among rotational capillary waves, an evolution equa- 
tion corresponding to the rotational case should be 
derived. The presence of non-zero vorticity invalidates 
the existence of a velocity potential for the flow, as is 
the case for irrotational flows, and harmonic function 
theory is not readily available for the analysis. Two 
common methods of analysis to find the exact form 
of the dynamic equations arc the method of multiple 
scales and variational techniques [10]. Both methods are 
applicable for flows with vorticity [3, 5, 9, 24]; this is 
work in progress. 
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